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Abstract. The low energy region of certain transition metal compounds reveals dramatic correlation effects 
between electrons, which can be studied by photoelectron spectroscopy. Theoretical investigations are often 
based on multi-orbital impurity models, which reveal modified versions of the Kondo effect. We present a 
systematic study of a multi-orbital Anderson-like model, based on a new semi-analytical impurity solver 
which goes beyond simple modifications of the well known NCA. We discuss one-particle excitation spectra 
and in particular the role of level positions and Coulomb-matrix elements. It is shown that the low-energy 
region as well as the overall features of spectra critically depend on the model parameters and on the quality 
of the approximations used. Recent photoelectron experiments and corresponding existing calculations are 
put into perspective. An interesting crossover scenario between different regimes of ground states with 
characteristically different local correlations is uncovered. 

PACS. 71.10.-w Theories and models of many-electron systems - 71.20.-b Electron density of states and 
band structure of crystalline solids - 71.27.-|-a Strongly correlated electron systems; heavy fermions - 
71. 55.-1 Impurity and defect levels 



> 

m 

00 

d 

OO 

o 



1 Introduction 

The Anderson model with a spin-degenerate local level hy- 
bridized to a broad conduction band has served in studies 
of the Kondo effect [1-3] and, combined with approxima- 
tions like DMFT, also in investigations of the behavior of 
lattice-periodic systems [4-6]. In connection with lattice 
systems it plays the role of an effective impurity which, 
as part of a self consisting cycle, has to be solved for the 
Greens functions of local electrons by use of appropriate 
impurity solvers. Among these impurity solvers are nu- 
merical ones like QMC [7,8] or NRG [9,10] and so called 
semi-analytical ones, which are based on direct perturba- 
tion theory with respect to hybridization and which rely 
on numerical procedures to solve implicit equations for 
propagators formulated via skeleton-diagram techniques 
[11-15]. The former essentially exact methods in fact can 
properly describe the low energy region of the model but 
are either limited to relatively small interaction values or 
high temperatures due to numerical obstacles and exhibit 
deficiencies for higher excitation energies. The newly de- 
veloped semi analytical methods like ENCA, SUNCA, and 
CAl [16-18], which improve the well known NCA [13,19, 
20, 14] via inclusion of vertex corrections with crossing di- 
agrams, on the other hand are well suited for the whole 
range of excitation energies. 

Capturing electron correlations due to strong local in- 
teractions in realistic transition metal compounds, in the 
framework of band structure theory or in form of quasi- 



particle band structures derived from Greens functions, is 
an active field of research [21-23]. Difficulties arise from 
the effective local building blocks needed in the theories, 
which in most cases are more complicated than a simple 
Anderson impurity. 

A connected problem of equal relevance arises in the 
interpretation of data from photoelectron spectroscopy 
(PES) in Ce-compounds, which hopefully can demonstrate 
the Kondo physics via direct observation of the formation 
of a many-body resonance near the Fermi level [24,25]. 
Here it is the influence of crystal-field levels and spin-orbit 
splitting, possibly also of higher levels on Ce with double 
occupancy, which complicates the picture and greatly en- 
hances the necessary numerical effort. Since the resolution 
of inverse PES still is not sufficiently high enough to mea- 
sure directly all details needed of the DOS above the Fermi 
level ep, one relies on direct PES together with extrapo- 
lations into a limited region of thermally occupied states 
above ep ■ Although a consistent picture of the many-body 
effects in some of the Ce-compounds arises [26] , some un- 
certainties remain. In particular, comprehensive and reli- 
able information about the whole spectral region of en- 
ergies up to the excitation energies of local states with 
differing occupancy and their influence on the low-energy 
region is highly desirable. 

Attempts of calculating one-particle excitation spec- 
tra for multi-orbital atoms in a metallic host in the frame 
of direct perturbation theory have so far been based on 
the NCA as a theory designed for the case of infinite lo- 
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cal Coulomb repulsion U [26,27], with some simple adap- 
tations for the additional orbitals [28,29]. Vertex correc- 
tions have thus been neglected and only two possibilities 
for the occupancy of the local shell arc realized. Published 
attempts to improve on this kind of approach have not fur- 
nished e.g. improved excitation spectra yet [30]. Neverthe- 
less, it is known, that even for the simple Anderson model 
important modifications arise from better approximations 
and more realistic values of the local Coulomb repulsion, 
concerning e.g. the characteristic low energy scales as well 
as the bulk of the more remote spectrum above the Fermi 
level [16]. 

A recent study [18] has demonstrated that among a 
few existing qualitative improvements of NCA, some of 
them with even somewhat higher acciiracy, the ENCA can 
serve as a rather reliable approximation in a large region of 
excitation energies, and that the numerical effort is limited 
and can be handled well. ENCA goes considerably beyond 
NCA; it includes all vertex corrections up to second order 
in the hybridization strength and describes the region of 
all finite U. ENCA e.g. correctly reproduces the low energy 
scale Tk with its dependence on model parameters [16]. 

We have generalized the ENCA to the case of a multi- 
orbital model and present here a systematic study of this 
new ENCACF. We will discuss some interesting features 
arising from variations of shell energies and different lo- 
cal Coulomb-matrix elements. As it will turn out, consis- 
tent interpretations of many features of the spectra cal- 
culated are possible, r(n'(>aling somc^ interesting physical 
trends and hopefully enabling an improved understanding 
of measured spectra. 

In the following section 2 we shortly present the model 
and the construction of our ENCACF approximation and 
we also compare some basic features with earlier approxi- 
mations. Section 3 contains our main results in form of cal- 
culated series of spectra for varying model parameters and 
relevant physical interpretations. A summary and some 
remarks on future developments concludes our paper in 
section 4. 



2 A multi-orbital Anderson model and the 
ENCACF approximation scheme 



When setting up realistic models for transition metal atoms 
with open inner shells, one is confronted with the problem 
of dealing with several free parameters which, though in 
principle to be determined consistently from first princi- 
ples, can usually not be fixed with sufficient accuracy. For 
definite host materials some information from experiments 
and from basic calculations like band structure theory is 
available and may help to define relevant regimes. In a 
situation, where physics is intricate and the validity of 
approximations is questionable, an approach via a some- 
what restricted model seems appropriate. We therefore 
will study the following generalization of the well known 



Anderson model: 



+ 



(1) 



+ + X! H ^J'^^l'^' ■ 

j j<l cr,a' 

Here ct^^.^ creates an electron in band j with crystal mo- 
mentum k, spin a and energy cjk, f]„ creates an electron in 
a local level j with spin a and energy Cj, and hja = fj^fja, 

fi'jkir = Cjf.^Cjk(7 are the corresponding occupation- number 
operators. The Vj denote (real) matrix elements for local 
hybridization between an electron in band j and one in 
a local level j, both with the same spin. Thereby we as- 
sume, as usual, that one-particle states can be classified 
according to their crystal symmetry with respect to the 
position of the impurity. With j being the relevant quan- 
tum number, hybridization in our model thus respects 
crystal and spin symmetry separately. Insofar as only the 
local hybridization intensities Aj{u}) = V^p'^^^uj), with 

(w) = J2k ^('^ ~ ^ik) being the DOS of band j, en- 
ter the following calculation, one may likewise speak of 
one single band and its respective symmetry components 
around the impurity. 

From a formal point of view, as well as from the physics 
of the model, conservation of j is important for the proper 
implementation of our diagrammatic approach, i.e. in the 
identification of irreducible pieces. We comment in the 
conclusion on how this is to be changed when studying lat- 
tice systems or complexes with several centers. In general 
the local Coulomb matrix elements Uji differ, depending 
on the states j and I involved. 

Altogether the model (1) contains a whole set of pa- 
rameters, including the functions p^^\lu). Unless one is 
able to derive all of these to a good accuracy with a spe- 
cific situation in mind sonic reasonable assumptions have 
to be made in order to restrict the wealth of detailed infor- 
mation obtainable to some typically situations and trends. 
We therefore assume just two local levels ei < £2, i.e. 
j = 1,2, without further degeneracies besides spin, and 
also set the two band dispersions eifc = e2k equal, as well as 
the corresponding hybridization matrix elements Vi = V2 ■ 
These are not completely unphysical assumptions; they 
can be realized in a highly symmetric situation, e.g for 
p-states aligned along different axes of a cubic crystal. 

For the transition metals we have in mind, the Uji will 
typical be of the order of the bandwidth W, and U12 might 
be somewhat smaller than Un and 1/22- We will consider 
the cases Uu = U22 = U12 = U and U\\ = U22 = U with 
varying f/12. 

Assuming a simple-cubic structure and a half-filled 
tight-binding band centered around the chemical poten- 
tial at efe = = 0, the unperturbed conduction band 
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DOS p^'^'^luj) stretches from uj = —3 to a; = 3 with van- 
Hove kinks at co = ±1. We work with a fixed Anderson 
width A = 7rF2p(o)(0) = 0.3 = W/20. Typically ei = -1 
lies in the occupied band region, Uji exceeds |ei|, and £2 
varies between ei and the upper band edge. Thus, a typ- 
ical Kondo regime in the sense of the original Anderson 
model is studied. 

The model (1) allows for a series of local occupation 
numbers, which in increasing order cause total shell en- 
ergies elevated by one of the Coulomb energies Uji. We 
restrict ourselves mostly to the three lowest values njoc = 
0,1,2 by properly truncating the system of integral equa- 
tions to be introduced below. This implies the following 
local shell-sates with the corresponding energies: 



n-ioc 


partial occupations 


energy 





nia- = n2a' = 


Eo = 


1 


nic = 1, ni-a- = n2a' = 


Ei-i = €1 


1 


n2a = 1, n2-<T = nia^ = 


El:2 = £2 


2 


nia = ni-a = 1- n2a' = 


£^2:11 = 261 + r/11 


2 


>l2a = Il2-a = 1. "la' = 


£-2;22 = 2(.2 + U22 


2 


n-la = n2a' = 1 

ni-cr = n2-a' = 


E2-12 = ei + £2 + Ui2- 



The unperturbed one-particle excitation energies involve 
all possible differences between two of the above levels 
with difference one in the total local occupation number 
nioc- This simple fact is useful to remember when inter- 
preting calculated spectra in the next section. For a com- 
plete classification of the high energy region the scheme 
has to be extended to nioc = 3 and n;oc = 4. This should 
be kept in mind, when spectral weight has to be accounted 
for, which is missing from the calculation due to the ne- 
glect of the corresponding propagators. 

The theoretical setup of semi analytical impurity solvers 
is based on the early work of Keiter and Kimball [11,31] 
and was brought into a form used today about 25 years 
ago [13,14,20]. It starts as "direct" time ordered pertur- 
bation theory in one-particle matrix elements between lo- 
cal and band states or in transfers between neighboring 
sites; Wick's theorem and Feynman diagrammatics can- 
not be used due to the presence of local interactions in the 
unperturbed Hamiltonian [12,32]. An alternative though 
equivalent approach uses additional, unphysical particles 
("slave bosons") in order to stay in the frame of Feynman 
perturbation theory [33,34]. 

In eitlica' formulation skeleton-diagram techniques are 
used to perform the infinite summations, which are neces- 
sary to capture the essence of the many-body effects, i.e. 
the Kondo effect as a prototypical example, determining 
the ground state and low-lying excitations. The implicit 
equations for various types of propagators of local many 
body states, formulated via skeleton diagrams, constitute 
a system of singiilar complex integral equations, which 
apart from a few special simple cases have to be solved 
numerically. Since the perturbational setup and the par- 
ticular diagram techniques are well documented [32,19, 
13, 12, 15], we will restrict our presentation to the new fea- 
tures and approximations developed for the multi-orbital 
model of equation (1). Various approximation schemes and 
relevant numerical techniques have recently been investi- 



gated and compared, regarding their quality and validity, 
in applications to the Anderson model and to simple lat- 
tice models [18]. We will rely on the experience gathered 
there and will adopt the ENCA in a proper generalization 
to the multi-orbital situation. Direct perturbation theory 
involves the bonus of allowing a direct interpretation of 
diagrams in terms of physical processes [32,12]. Therefore 
we will present corresponding pictures in order to shortly 
explain the ideas behind the ENCACF used in the follow- 
ing calculations. 

Propagators for states of the local shell evolve via ex- 
citation processes along a time axis with imaginary time 
for technical reasons. Pieces of processes, which are inter- 
nally connected via lines visualizing excitations into the 
band, constitute irreducible self-energies for the propaga- 
tors. The self-energy diagrams corresponding to a local 
state with just one electron of spin a in orbital £1 are de- 
picted in Fig. 1(a). The earliest, i.e. lowest band electron 
line, labeled (fc, u'), determines via its direction, upward 
or downward, whether the following local state has zero 
or two electrons. In the first diagrammatic contribution 
it is the vacuum state |0), in the following two cases the 
incoming electron, the spin of which is preserved in the 
hybridization event according to equation (1), is absorbed 
either in level €1 or in level £2. Two electrons in level £1 
are to have opposite spins, i.e. a' = —a here, while in 
the last contribution the second spin-direction a' is free 
to be summed over as well as the crystal momentum k. 
All later hybridization events can be collected in the ver- 
tex functions, represented by the different triangles for 
the three cases. Approximations to be made concern these 
vertex functions. In a simple non-crossing approximation, 
adopted to this multi-orbital case (SNCACF), the trian- 
gles are just replaced by bare vertices. 

Although this scheme looks very simple, it is the use 
of full propagators for all local linc!s shown, with the self- 
energies under debate inserted, which already makes it 
highly non-trivial and allows e.g. for a qualitative descrip- 
tion of Kondo physics. Apparently, in a picture with prop- 
agators for the bare states, series of excitations enclosing 
each other without crossing are generated. The ENCA- 
scheme introduces the vertex corrections up to second or- 
der in the hybridization and thus incorporates many more 
processes with crossing band-electron lines. For the An- 
derson model it was shown to improve considerably on 
the SNCA, reproducing e.g. the correct Kondo scale [16]. 

In the generalization to our multi-orbital model the 
ENCACF incorporates vertices like the one shown in Fig. 1(b). 

This example suffices to clarify the building principle 
of the approximation. The other two vertex functions ap- 
pearing in Fig. 1(a) are constructed in an analogous fash- 
ion. They depend on two energy variables, one of them 
complex, and have to be treated in all analytical equa- 
tions with great care; numerically, their inclusion consid- 
erably enlarges the calculational effort. We will not make 
explicit all the diagrammatic self-energy equations needed 
and the corresponding additional set of diagrams neces- 
sary to build up the one-particle- Green function consis- 
tently in ENCACF. In principle, the diagrams are ob- 
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Fig. 1. Diagrammatic representation of ionic self-energies (a) and vertex corrections (b) for a single-impurity Anderson model 
(SIAM) in direct perturbation theory. Physical processes arc arranged vertically along an imaginary time axis (broken line) 
which bears an energy variable z after Laplace-transformation. Presence of an electron in the local shell is indicated via a wiggly 
line on this time axis. Excitations of band electrons (straight lines) take place at hybridization vertices (dots on the time axis). 
All line crossings in these diagrams can be attributed as a correction to one of the vertices, which axe drawn as triangles. 



tained from the ones of Fig. 1 and the additional ones 
not shown, by cutting one band electron line. Following 
the remarks above and the well known recipes for a tran- 
scription of such diagrams into formulas [32,12] it is a 
straightforward task to complete the systems of integral 
equations for local propagators and for defect propagators 
[35] , as well as to assemble all contributions to the Green 
functions and the excitation spectra derived from them. 

We have solved the full ENCACF equations to con- 
vergence and will present one-particle excitation spectra 
derived from Green functions according to 



(i,^ = l,2) 



(2) 



In most cases we will present the complete spectrum 

j 

which in our case is normalized to 2, since it contains ex- 
citations of both local levels ei and £2 with definite spin a. 

First we compare a typical calculation with some simpler 
cases, which in part are known from the literature. 



SNCACF with the same value oi U — 2 reveals this still 
to be true in the multi-orbital model: In fact, the spectral 
weight accumulated near the chemical potential at a; = 
is much larger in the first case, although the two different 
temperatures likewise approximately reflect the respective 
approximations to the Kondo scales of the model without 
the second orbital (see inset of Fig. 2). 

A comparison with the SNCACF-spectrum at larger 
[/ = 4 demonstrates the shortcomings of a simplified cal- 
culation, where all doubly occupied states are projected 
out via letting U go to infinity: Considerable spectral weight 
is driven out of the band region, in form of increasingly 
narrow spikes. The ENCA-calculation without the second 
orbital, finally, helps to identify the additional spectral 
weight induced by the level 62, not only in the region of 
high excited states around w « 2, but also near oj = 0, 
where the many-body resonance acquires two satellites 
above and below. In a qualitative sense this behavior is 
known from earlier calculations using SNCA and U = 00 
[36,37]. Although the number of one-particle transitions 
contributing to these spectra is rather small due to the 
uniform value of Coulomb-matrix elements and to the re- 
striction to the subspace with rijoc ^ 2, is it worthwhile 



In Fig. 2 we present the a ENCACF-spectrum, a SNCACFfeo identify the bare values of transition energies. These 



calculation with the same and an elevated value of the 
local repulsion, together with an ENCA-spectrum of the 
simple Anderson model without the second local level £2 = 
—0.5. The Coulomb-matrix elements Uu = U22 = U12 
all have the same value U = 2, except for the second 
SNCACF -case where C/ = 4. The inverse temperature is 
/? = 50 for the ENCA-cases and /? = 200 for the SNCACF, 
to be comparable, since SNCA underestimates the many- 
body scale considerably. Comparison of the ENCACF and 



are for U = 2: Zie^°^ = 

Ei,2-Ea = €2 = -0.5, 



Ae 



(0) 



Ei-i 







ei 



-1, Aei"^ = 



Eo 



-Ei,i = e2 + U 



— E2:ll — El-l — €l - 

= 1.5, Aef = E2;i2 



f/=l, 

-Ei-2 = 



ei [/ = 1 and Ae)^> = E2-22 - E1.2 = €2 + U = 1.5. Ap- 
parently in a region of energies uj between +1 and +2 con- 
siderable spectral weight is accumulated. The resonances 
acquire widths on the one-particle scale A and in excess 
thereof via the local interactions and can therefore in part 
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Fig. 2. Local onc-particlo excitation spectrum of the two-orbital model (first three curves) and of a corresponding Anderson 
model with one orbital, calculated within SNCACF and ENCACF, respectively. The temperature T = = 200kg 

the two SNCACF calculations is chosen higher than in the two ENCACF cases, so that it becomes comparable with the 
approximate Kondo scale in each case. Other model parameters are ei = —1, €2 = 0.5, A = np^°\0) = 0.3 (Anderson width) 
and ?7ii = U22 = U12 = 3; energies are in units of the 3d=6*'' part of the width of a simple cubic bandstructure. 



not be well resolved individually. The resonance originat- 
ing from Aef^ , i.e. the lowest peak in the figure is known 
to be broadened by 2 • z4 in a pure ENC A-calculation with- 
out the higher local level due to the blocking of the res- 
onance, if an electron with opposite spin already resides 
there [36] . This blocking effect apparently is enhanced by 
the presence of the second local level. 

The visibility of excitation peaks in these spectra is 
largely determined by the thermal occupation probabil- 
ity of the initial and final shell state. Whereas e.g. the 
transition Ae^^^ possesses a weight (Pq) + (A;1ct) ~ 5 
((Pjvf) being the thermal occupation number of the local 
shell-state M with Pm the corresponding projection op- 
erator), the next transition at ^€2'^ has negligible weight 
(Pq) + {Pi:2a) 1 in the regime considered here, with sin- 
gle local occupation in level e\ being by far most probable. 
The high- lying peaks in the U = 4-case shown in 2, like- 
wise have weight near 1, but already indicate a crossover 
to narrow spikes, since broadening begins to decrease due 
to a lack of nearby band states. In this way, the rough fea- 
tures of single particle excitation spectra can completely 
be classified. 



3 Discussion of one-particle excitation 
spectra and correlation effects 

Figs. 3 and 4 are designed in order to identify correlation 
effects in the spectra, and in particular the influence of 
the presence of the second local one-particle level with 



higher energy €2 > ei . A clear hint to the important role of 
many-body correlations, beyond the interaction-induced 
splitting of the local levels into a lower part at about and 
a higher part near tj -\- U both with fractional weights, is a 
pronounced temperature dependence in a certain interval 
of energies w. 

Fig. 3 shows sequences of spectra for varying temper- 
ature, typically reaching from about a third of a charac- 
teristic low energy scale fc^T* (/3 = 100) up to about ten 
times ksT*. A value of about fcsT* fa can be deduced 
from the half- width of the narrow peak at w = 0; it would 
coincide with the Kondo scale fesTx, if the higher local 
level were not involved. Nevertheless, the choice of param- 
eters for the figures implies a situation like the one found 
in the classical Kondo regime, with €2 = —0.8 or —0.5 
clearly larger than ei = —1, both negative, but Un = 
U22 = U12 large enough to favour a ground state occupa- 
tion nioc = X)j=i 2 S<7 ~ both cases of Fig. 3, for 
£2— ei = 0.2 small (part a) and for e2— ei = 0.5 larger (part 
b), the curves in a narrow region around the chemical po- 
tential Lo = 11 — prove to be strongly temperature depen- 
dent, whereas the remote parts of the spectra are much less 
affected by temperature variations. The qualitative pic- 
ture, known from earlier work, of a three peak many-body 
resonance forming when T* is approached from above is 
clearly seen in Fig. 3(b). In Fig. 3(a), however, the low- 
est and weakest of these peaks is concealed as a very 
flat shoulder on the left side of the main resonance in 
the middle. This result again demonstrates the stronger 
lifetime effects seen in the ENCACF as compared to a 
simpler SNCACF-calculation, which tend to broaden the 
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Fig. 3. Temperature dependence of the local one-particle excitation spectrum of the two-orbital model in ENCACF- 
approximation. Parts (a) and (b) differ in the value of the second local particle level €2 = —0.8 and 62 = —0.5, respectively. 
Other model parameters are as in Fig. 2. 



resonances further and increase the scale Tk considerably. 

The origin of the thrcc-pcak structure has been explained 
consistently [36] : The ASR of the Kondo scenario acquires 
two satellites, one on each side, due to additional local ex- 
citations in the spectral region shaped by the many body 
effects. A local excitation from level ei to level £2 pre- 
serves the possibility of forming a Kondo singlet of dif- 
ferent kind, and the corresponding new ASR is shifted 
via the foregoing excitation process by Ae = €2 — ei above 
the chemical potential. The inverse process leaves its trace 
near oj = —Ae; its thermal weight however, is reduced by 
the small occupation probability e2 compared to the lower 
level €1. 

Fig. 4 demonstrates more clearly that this interpreta- 
tion in fact holds true: The two side peaks of the main 
resonance actually move with the position of the higher 



local level £2 and correspondingly also become less pro- 
nounced with growing distance. The higher energy satel- 
lite can nevertheless be followed through the high energy 
region around and above ei + U = 2; it appears in the 
curve for £2 = 2 on the right hand side of this peak, but 
below the peak at energy 62 + U = 5. The choice of pa- 
rameters for this case according to the listing of excitation 
energies in the last section, only give unperturbed peak 
positions at Ae^°^ = ^ = -1, Aef'> = Ae^°^ = Z\£^°^ = 

€2 = ei + U = 2, Zi£f = Aef = e2 = €2 + U = 5, so 
that the given identification of the feature near w = 3.3 
seems to be unique. The temperature chosen via /? = 50, 
for Fig. 3 and all the following ones corresponds to T* 
or somewhat lower, so that all many-body effects should 
clearly be visible. 
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Fig. 4. Dependence of the local one-particle excitation spectrum of the two-orbital model on the position of the second level 
£2 in ENCACF-approximation. The temperature compares roughly with the Kondo scale, other model parameters are as 
in Fig. 2. 



Also the energy region below the chemical potentials 

reveals some interesting features. In particular, the way in 
which the exact form of the lower main resonance of the 
original Anderson model without the higher local level 
is reconstructed with growing £2, allows to identify the 
amount of additional scattering due to the blocking be- 
tween electrons in levels ei and £2. Even hero, the lower 
and weaker side peak can be traced for a limited region of 
increasing e2-values as a shallow flank, although its ther- 
mal weight decreases and it is obscured somewhat by the 
changing form and position of the main resonance aris- 
ing from Aef'^ = ei = —1. The increase of the center of 
this broad main resonance with growing £2 reflects the loss 
in ground state energy , when level interactions become 
weaker. The shrinking width of this resonance indicates 
the loss of blocking of level ei due to the presence of elec- 
trons in level £2, when this energy increases. 

The next figures are devoted to a study of the influ- 
ence, which the magnitude of local Coulomb-matrix ele- 
ments has on the spectra and on the correlations. In Fig. 
5 we keep the equality Un = U22 = U12 ^ U, but vary 
U in integer steps from 2 to 5. Unlike the variation of 
temperature in Fig. 3 or the variation of the upper level 
position in Fig. 4 the variation of U does affect the width 
of the central many-body resonance. This T* gets smaller 
with increasing U in accord with the behavior of Tk in the 
Kondo effect. The decrease of T* is moderate and effec- 
tively stops when the upper resonances ei + U and e2 + U 
cross the upper band edge at lu = 3. The splitting of these 
two highest peaks is flxed approximately at our value of 
£2 — £1 = 0.5, with £1 = —1. The width of both peaks 
decreases from U = 3, in which case they just stay inside 
the band region, to U = 5, where both are situated above 
it. Increased scattering inside the region with high density 
of band states washes these peaks out which best can be 



seen in the region 1 < w < 2.5 for the lowest value U = 2. 

The height of the central many-body resonance changes 
only slightly, showing that we stay in the near saturated 
region T<T*(C/). 

A very interesting perspective opens when variations 
of the Coulomb matrix element U12, i.e. the repulsion be- 
tween electrons in different local levels £1 and £2, are per- 
formed keeping Un = U22 fixed. Fig. 6 shows the cor- 
responding spectra, again for fixed £1 = —1, £2 = —0.5, 
(3 = 50 and Uu = U22 = 3, and for a variation of U12 
in integer steps between and 5. It is apparent, that a 
considerable transfer of spectral weight from above to be- 
low the chemical potential takes place when going from 
values U12 > 1 to values U12 ~ 0. In fact, the condi- 
tion £2 + U12 = 0, i.e. U12 = 0.5 marks a crossover here: 
For U12 > 0.5 a ground state with single level occupation 
nioc = 1 is stable in zeroth order as in all cases before, 
whereas for U12 < 0.5 the stability goes over to a new 
ground state with nia = n2a = \, i.e. nioc = J2ja ^jc = 2- 
This changes drastically the thermal occupation probabil- 
ities of local shell states and the corresponding transitions 
as well as the nature of the many-body correlations. 

Generally, the thermal weight of a one-particle tran- 
sition from shcU state M to M' is (Pm) + (Pm'), i-c 
the sum of the occupation probabilities if the initial and 
the final state, in zeroth order. The corresponding exci- 
tation energy is Ae = Em' — Em and n^^^ = nioc + 1 
by definition. For large U12 and £1 < £2 in the regime 
considered here, (A;1<t) is approximately 5 and all other 
{Pm) are small. Consequently, the transitions with ener- 
gies Aef^ = £;i;i-Po = ei,^er = £^2:ii-Pi;i = £i + C/ii 
and A<^^^^ = i?2;i2 — Ei-i = £2 + U12 are strongly and 
the others are wciakly represented in the spectra. For i.e. 
f/i2 = 4 this explains the 3-peak structure in the high 



8 



Grewe et. al: Multi-orbital Anderson models and the Kondo effect 



1 



0.5 



1 1 1 

U^2. — 


1 




1 ii 


1 ;l 


1 




U = 3. 










1 

1 1 




U = 4. — 






1 i i 




h 
i i 




- U = 5. — 










1 1 
1 

1 1 
1 1 








/A 1 
/ \ ' * 1 
^ \ ' \ ' 


1 ' ' 




1 1 
1 

1 1 
1 1 
1 1 

1 1 




J 




\ / 
\ / 

A 

/ \ 

y \ 


\i' \ / 




; > 1 
; ' ' 

', ' 1 
'\' \ 


I ; 
1 \ 
1 ^ 
1 V 

\ 




1 






-■,1-.- 


-—4-—'::: 



Fig. 5. Dependence of the local one-particle excitation spectrum of the two-orbital model on the Coulomb matrix element U 
for the case Un — U22 = U12 in ENCACF-approximation. The temperature and second level position are and £2 = —0.5, 

respectively. Other model parameters are as in Fig. 2. 




Fig. 6. Dependence of the local one-partiele excitation spectrum of the two-orbital model on the particular size of U12 of the 
Coulomb-repulsion between an electron in level ei and one in level 62 in ENCACF-approximation. Other model parameters are 
as in Fig. 2. 



energy region with peaks shifted somewhat outward from 
Ae^°^ = -1, Aef^ = 2, and Ae^^ = 3.5. For small U12, 
on the other hand , {P2:i<72a'} is approximately j and the 
other probabilities arc small, so that only the transitions 
with energy Aef^ and /Acg"^ = £2-12 — Ei-2 are strong 
and the others weak. The values i.e. for U12 = 0.5 are 
Ae'f^ = and Ae^^ = —0.5, for U12 = one obtains 
Ae'f^ = -0.5 and Ae^^'^ = -1. 

These considerations ah'cady give a rough explanation 
for the observed shift of spectral weight, if one takes into 



account that peak positions are usually rcnormalizcd away 
from the chemical potential by interactions. The crossover 
region itself bears interesting physics in it. At the point 
U12 = 0.5, ei = —1 and £2 = —0.5 as before, the two states 
Ei-i = ei = —1 and £2-12 ~ ei+e2 + C^i2 = —1 with single 
and double occupancy arc degenerate ground states in zc- 
roth order, which puts the system in some kind of interme- 
diate valence regime. This corresponds well with a broad 
distribution of spectral weight near lu = 0. The effect is 
seen in the curve for U12 = 1, and even more pronounced, 
in the following Fig. 7 for U12 = 1, with ei = e2 = — 1. 
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Fig. 7. Like Fig. 6, but for ei = €2 = —1, i.e. the two local levels ei and €2 are equally occupied and contribute symmetrically, 
i.e. to separate Kondo effects at U12 = 0. 



Consistently, in Fig. 7, for vanishing J7i2 = a marked 
Kondo effect with a thin ASR near w = reappears. In 

this case the ground state -E2:i2 = —2 is well separated 
from the first excited states Ei-i ~ i?2;2 — — — 1- 

With respect to the local many-body correlations the 
crossover scenario implies the following consequences: f/12 
can be viewed as a measure for the competition between 
tendencies to occupy each of the two local levels. For 
large C/12 > Un this competition favors the lower local 
level ei and its corresponding Kondo effect. Level €2 then 
is only populated with a small probability and furnishes 
only small side peaks near the ASR. Correlations due to 
a virtual Kondo effect of level £2 stay incomplete. With 
U12 ~ Uii = U22 reaching the value of the other Coulomb- 
matrix elements all of the different forms of local double 
occupancy become energetically comparable, so that the 
side peaks are strengthened somewhat due to increased in- 
terference of the correlations. Below the crossover to the 
doubly occupied ground state, however, individual Kondo 
effects for level ei and £2 start to emerge. They do not show 
up for small absolute values of ei and 62, where ground 
states with different local occupancy have comparable en- 
ergies (on the scale of the Anderson width) and low- lying 
charge fluctuations cause a broad peak near a; = 0. If 
ei and 62 are both well below the chemical potential a 
complete decoupling then occurs at the point U12 = 0, 
where each of the levels only interacts with its own band 
of the same symmetry. Due to different level positions the 
two Kondo temperatures Tkj ^ cxp{ — y—^jy^) , Jj « 

+ e +u come out largely different, as a pos- 

sibility. With e.g. Tk^ <T < Tk^ the ASR due to €2 is 
completely formed, whereas the one due to ei is still in 
an early stage; the low-energy excitations then would be 
largely determined by the Kondo effect of the higher lying 
local level. 



For our set of parameters Tk^ should be considerably 
higher than Tk^, and the position of £2 = —0.5 indicates 
an intermediate valence region. The corresponding "ASR" 
then significantly broadens and acquires some degree of 
linewidth via charge fluctuations and simple one-particle 
scattering oc V"^. Thus due to the role of temperature 
variations relative to Tr-j and Tk^ and due to different 
level positions €2 in the interval between ei and the chem- 
ical potential /x = 0, the crossover scenario presented here 
can involve quite different shapes of the one-particle ex- 
citation spectrum. In particular, when even more compli- 
cated multi-orbital models are used e.g. in calculations of 
correlated bandstructures, great care seems necessary to 
understand and control all parts of the calculation and 
the physical relevance of all features. We will elaborate 
on this cautionary remark below also in connection with 
some technical aspects of the approximations used. 

With Fig. 8 we present further evidence to support 

our picture of the crossover scenario. We have resolved 
spectra from Fig. 6 on each site of the crossover, into the 
partial spectra, which incorporate only the excitation on 
one of the local levels. Fig. 8 clearly demonstrates the 
strong dominance of the Kondo effect on the lower level ei 
when U12 is large: The ASR seen for Ux2 = 4 is completely 
contained in the partial spectrum p-i^{Lo), whereas P2cr(w) 
is continuously flat near w = 0, although T = -^^^ is of 
the order of a hypothetical Tr-^ of isolated level £2- 
In Fig. 7 on the other hand, we have seen a symmetrical 
situation with ei = £2 = — 1. The ASR, which shows up 
for U12 =0, is clearly shared in equal parts by both partial 
spectra. A similar effect in the spectrum with lowest C/12 of 
Fig. 6 is not seen due to the different nature of correlations 
in the Intermediate Valence regime. 

The regime with the doubly occupied local ground 
state is apparently not as well described by our ENCACF- 
approximation (not to mention the simpler SNCACF) as 
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Fig. 8. Partial local one-particle excitation spectra , P2ct(w) in comparison with the combined spectrum pa{oj) = 

PictC^) + P2ct(w). The model parameters are as in Fig. 6 with U\2 = 4. 




Fig. 9. Local one-particle excitation with restricted (local occupancy < 2) and full (local occupancy < 4). The model parameters 
are as in Fig. 7 with U\2 = 0. 



the former regime with the singly occupied local ground 
state. A plain hint comes from the lack of spectral inten- 
sity in the high energy region above w = 0. More subtle is 
a corresponding influence on the correlations determining 
the spectrum nearby the chemical potential. With a local 
ground state derived from £^2:12 = ei + £2 + U12 and U12 
small, excitations to a local state with triple occupancy 
and energy ei + 2e2 + U\_2 + U22, i-e. Ae^^ = £2 + U22, or 
energy 2ei+e2 + C/i2 + f/ii: i.e. Ae^^^ = ei + Un, contribute 
to the dynamics. Their importance compares with that of 
Aef'^ = 62 + U12 and Ae'^^ = ei + t^i2, similar to the im- 
portance of both one-particle excitations at ei and ei + U 
for an Anderson model with finite U. This deficiency can 



only be cured, when the Hilbert space underlying the cal- 
culations is enlarged to account for triple ( and quadruple) 
local oc:c;upancy. Thereby the system of integral equations 
for propagators again becomes considerably enlarged. 

We have also solved the correspondingly completed 
ENCACF-theory numerically and present a spectrum in 
Fig. 9. It is calculated for the most interesting case of 
ei = £2 = -1, U12 = and Un = U22 = p = 50 and 
is compared with the incomplete spectrum discussed be- 
fore. The difference, particularly in the high-energy region 
is clearly visible. On the othc;r hand these differences for 
cases, where the local ground state is singly occupied, es- 
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sentially only show up at very high excitations energies 
u) > 2U and can safely be ignored. 

A second technical remark also refers to our note of 
caution made above. Although from a puristic point of 
view and in agreement with physical insight, the system 
for Ui2 = decouples dynamically into two independent 
ones, this is in general not realized in approximate version 
of direct perturbation theory, which from the beginning 
incorporates C/12 into the unperturbed local system. For 
example, with U12 = and ej < —Aj < < Aj < Cj + Ujj 
{j = 1,2), one should observe two independent Kondo 
effects. Although as wc have demonstrated, this is actu- 
ally found for large enough absolute values of both ej , we 
are not so confident about our findings in the crossover 
regime. The point can be visualized with Fig. 1(a). The 
third self-energy diagram shown there has an intermedi- 
ate state ni = 712 = 1, which is the ground state for small 
C/12. Fluctuations like the one shown via the downrunning 
band electron line, as projected on the dynamics of level 
ei, can only take place when an electron is present in level 
€2- This establishes a spurious kind of correlation, which 
only is eliminated from the calculation when all of the pos- 
sible vertex corrections are actually taken into account in 
the vertex shown at the top. Only this can lead to all com- 
binations of individual time orders, which together would 
decouple the two subsystems. 

Our ENCACF captures the leading part thereof via 
the vertex corrections shown in figure 1(b). It is, however, 
not settled yet how important higher contributions with 
more crossings are in the different regimes. Although there 
is some experience for the simple Anderson model [38,37, 
18], suggesting a qualitatively well behaved ENCA, in par- 
ticular the crossover regime found here might need even 
better approximations. 

4 Conclusion and Perspective 

Our investigation of a multi-orbital model with a new 
crossing approximation, the ENCACF, which goes beyond 

a simple modification of NCA and works for all finite val- 
ues of Coulomb parameters, has revealed a variety of pos- 
sible forms of one particle excitation spectra. We have 
studied the role of different model parameters, tempera- 
ture T, energies ej of the two local one-particle levels with 
different symmetry, and of the Coulomb-matrix elements 
Uji between electrons in level Cj and e;. Different regimes 
could be characterized, according to the nature of the local 
ground state and of the local correlations. These regimes 
exhibit different modified forms of a Kondo effect, when 
the levels ei and 62 are below the chemical potential and 
representative values of Uji are large enough. 

Single particle transition between local states on the 
one hand and traces of many-body correlations on the 
other hand could both be identified in the calculated spec- 
tra and were discussed in a consistent framework. The 
crossover between regimes of different local groundstate 
occupancy involves interesting Intermediate Valence phe- 
nomena, where charge fiuctuations intervene. We have poin 
to the intricate physics, which determines the crossover 



region, and to some problems of approximations which 
possibly render a complete understanding difficult with 
the methods developed so far. Particular caution seems 
necessary when interpreting e.g. experimental data from 
photoelectron spectroscopy with insufficient approxima- 
tions, using e.g. simpler forms of direct perturbation the- 
ory like NCA, combined with the limit of infinite local 
Coulomb repulsion. Whereas we believe, that some cor- 
rect qualitative information can be gathered in this way 
and a correct overall picture be developed, our new cal- 
culations show, based on an essentially improved theory, 
show that quantitative fits may be very misleading: The 
relevant many-body scale, for example grossly differs and 
important spectral information from the high energy re- 
gion gets lost. Moreover, even the detailed spectral inten- 
sities in the low energy region depend on the rest of the 
spectrum and may come out incorrectly, in particular in 
real physical parameter regimes of transition metal ions 
which do not allow for an argument of universality. Nev- 
ertheless, our results also support some earlier findings, 
in particular for the spectroscopy of metallic compounds 
with Ce-atoms. For the use of multi-orbital models of the 
kind studied in connection with self- consistent methods 
for correlated band structures, on the other hand, our in- 
vestigation furnishes a rather poor perspective due to the 
complexity of the problem and the insufficient quality of 
suitable impurity solvers proposed so far. Here we hope 
that generalizations of our ENCACF may be helpful. 

In any case, multi-orbital impurity solvers for lattice 
systems with strong local correlations in most cases should 
be planned in an even more general form than our model 
of equation (1). The reason simply is that symmetry de- 
compositions of band states around different sites as cen- 
ters are not compatible : They combine the Bloch func- 
tions in a different way. As a consequence, the medium 
around a given representative site will mix the symmetry- 
decompositions, i.e. eife and e2k in our model. In that case 
the possible diagrammatic processes make up a larger set 
as before, with the consequence that the irreducible self- 
energies of ionic propagators become a matrix. This ma- 
trix involves also elements which are non-diagonal in the 
level indices j. The same complication arises in e.g. the 
case, where an impurity complex made up of two sites 
is inserted into a perfect crystal with orbitals on each of 
the two sites. This again spoils the conservation of a j 
quantum number of local electrons and leads to the same 
consequences as before. We have taken up studies of such 
generalized multi-orbital models and will report on results 
in a forthcoming paper. 
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